Parametrized tests of post- Newtonian theory using Advanced LIGO and Einstein Telescope 
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General relativity has very specific predictions for the gravitational waveforms from inspiralling compact 
binaries obtained using the post-Newtonian (PN) approximation. We investigate the extent to which the mea- 
surement of the PN coefficients, possible with the second generation gravitational-wave detectors such as the 
Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) and the third generation gravitational- 
wave detectors such as the Einstein Telescope (ET), could be used to test post-Newtonian theory and to 
put bounds on a subclass of parametrized-post-Einstein theories which differ from general relativity in a 
parametrized sense. We demonstrate this possibility by employing the best inspiralling waveform model for 
nonspinning compact binaries which is 3.5PN accurate in phase and 3PN in amplitude. Within the class of 
theories considered, Advanced LIGO can test the theory at 1.5PN and thus the leading tail term. Future obser- 
vations of stellar mass black hole binaries by ET can test the consistency between the various PN coefficients 
in the gravitational-wave phasing over the mass range of 1 1-44Mq. The choice of the lower frequency cut off 
is important for testing post-Newtonian theory using the ET. The bias in the test arising from the assumption of 
nonspinning binaries is indicated. 

PACS numbers: 04.30.Db, 04.25.Nx, ()4.8().Nn, 95.55.Ym 



I. INTRODUCTION 

General relativity is tested to unprecedented accuracies in 
the weak-field and strong-field regimes (see Ref. lit] for a 
review). From a theoretical perspective, tests of this nature 
were possible due to physically motivated but structurally 
simple parametrizations of the observable quantities which 
could have different values in different theories of gravity. In 
the weak-field regime, solar system bounds were largely as- 
sisted by the parametrized post-Newtonian (PPN) framework 
(see for example Ref. 12]). PPN formalism parametrizes the 
deviation of a general metric theory of gravity (with symmet- 
ric metric) from the Newtonian theory in the weak-field limit 
in terms of 10 free parameters up to order u/c, where v is the 
characteristic velocity of the object. General relativity (GR) 
is a special case of this class with specific values of these pa- 
rameters. 

Binary pulsar tests, which dealt with stronger gravitational 
fields involving compact objects but typical velocities of or- 
der V ~ lO^^c, probed GR in the strong-field radiative regime. 
The binary pulsar tests were performed effectively with the 
use of the parametrized post-Keplerian (PPK) llJ-la] formu- 
lation of the pulsar timing formula. The timing formula can 
be expressed as a function of Keplerian and post-Keplerian 
parameters, each one of which is a function of the compo- 
nent masses of the binary. A measurement of two of these pa- 
rameters enabled the determination of the individual masses. 
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The measurement of a third parameter would constitute a test 
of the theory, by requiring a consistency of the component 
masses in the mi - 1112 plane. Depending on the number of 
these PPK parameters that can be measured from the timing 
data of a binary pulsar, it enables many tests of GR (measur- 
ing n parameters allow « -2 tests). Binary pulsar observations 
also confirmed the quadrupole formula for the generation of 
gravitational waves. 



A. Gravitational waves and tests of GR 

The detection of gravitational waves (GWs) would be the 
first direct test of the consistency of gravitation with the prin- 
ciples of special relativity and would probe general relativity 
beyond the quadrupole formula ]|7|]. A subsequent detailed 
study of the properties of GW would next allow one to as- 
sess the validity of GR in the strong-field radiative regime. A 
prominent class of GW sources is compact binaries: neutron 
stars (NS) and/or black holes (BH) moving in circular orbit 
with velocities v ~ 0.2c. Within GR, using diff'erent analytical 
and numerical schemes, gravitational waveforms from these 
systems can be computed with very high accuracy 18]. Avail- 
ability of such high-accuracy waveforms will allow the appli- 
cation of matched filtering techniques to search for these sig- 
nals in the data from the GW interferometers such as LIGO ]|3] 
and Virgo ifioll. 

Beyond the detection of GWs, one would like to know 
whether one can perform tests of GR with the detected signals. 
Despite the use of GR waveforms in matched filtering (which 
essentially assumes that GR is the correct theory of gravity), 
several authors have argued that GW observations can be used 
to test GR and put bounds on various parameters in alternative 



theories of gravity. One of the first proposals towards testing 
nonhnear aspects of GR using GWs was due to Blanchet and 
Sathyaprakash, who discussed the possibility of measuring the 
"tail" eff'ect in the GW phasing formula IIIIII12I1 . Ryan pro- 
posed a method to measure various multipole moments of a 
binary system [1311 from the Laser Interferometer Space An- 
tenna (LISA) observations of extreme mass ratio inspirals. 
Will obtained the additional contributions to the GW phas- 
ing formula in Brans-Dicke theories Ill4ll and massive gravi- 
ton theories lUSll as a one-parameter deviation from GR and 
discussed the bounds possible on these corresponding param- 
eters from GW observations. These ideas were elaborated on 
in greater detail in a series of papers III6I - I21I1 . studying how 
various physical effects in the binary aff'ect the bounds. It is 
worth noting that these bounds possible on massive graviton 
theories will be complementary to those which are obtained 
by binary pulsar observations (see e.g. 12211 ). 

The basic idea of our proposal can be viewed as a general- 
ization of some of the existing proposals to test specific theo- 
ries of gravity like Brans-Dicke or massive graviton theories. 
General relativity and any of its parametrized variants have 
different predictions for the PN coefficients i/', in the phasing 
formula (for details see Sec. Ill CI ). Hence the accuracies with 
which the PN coefficients of GR can be measured translates 
into bounds on the values of these coefficients in any other 
theory. This leads to the question of how well can these coef- 
ficients be measured? 

One way to answer this question is to rephrase it as a param- 
eter estimation problem and measure each of the PN coeffi- 
cients, treating them all as independent of one another. Recall 
that, for nonspinning binaries, each one of them is a function 
of only the two component masses and hence just two param- 
eters are enough to describe the phasing formula up to any PN 
order Hence, if we want to treat each one of them (8 in all 
for the restricted waveform (RWF) at 3.5PN order) indepen- 
dently, there will be large correlations among the parameters. 
Our earlier work ll23ll has shown that this method works well 
only for a narrow range of masses for which the signal-to- 
noise ratio (SNR) is high enough to discriminate these terms. 
The masses are typically of a million solar mass and hence 
detectable by the LISA 111 HI. 

Since the high correlations among the PN parameters are 
responsible for the ineffectiveness of the above test, we ex- 
plored other possibilities which, though less generic, are vi- 
able, interesting and capture the essential features of the test. 
One possibility is to use a smaller set of parameters comprised 
of the PN coefficient to be tested together with any two of the 
remaining PN coefficients chosen as basic variables (to reex- 
press and parametrize the rest of the PN coefficients) ||25ll . 

While all these were based on inspiral waveforms, there are 
proposed tests based on merger and ringdown waveforms of 
the binary as well by measuring very accurately the various 
ringdown modes of the GW spectrum ll26l - l28ll . Keppel and 
Ajith 112911 revisited the bounds on massive graviton theories 
including the merger and ringdown contributions. Alexander 
et al pointed out that the GW observations can be used as a 
probe of effective quantum gravity which predict amplitude 
birefringence of the spacetime for the propagation of the GW 



signals yCj- Molina et al investigated the possible imprints 
of Chern-Simon theory of gravity in the GW ringdown sig- 
nals and its detectability with GW interferometers Bill . In 
brief, GW measurements can lead to interesting tests of vari- 
ous strong-field aspects of gravity. 

Recently, Yunes and Pretorius discussed a generalized 
framework called the parametrized post-Einsteinian (ppE) 
framework to describe various fundamental biases in theo- 
retical modeling of GW and express them in a parametrized 
manner Il32ll . They used the existing knowledge about vari- 
ous alternative theories of gravity such as Brans-Dicke, mas- 
sive graviton theories and Chern-Simon theory to write down 
a generic Fourier domain gravitational waveform, which is 
parametrized in terms of a set of amplitude and phase vari- 
ables. They also considered the contribution from merger and 
ringdown phases of the binary evolution beyond the inspiral. 
This parametrization in the inspiral regime can be considered 
as a generalization of our earlier proposal in Ref. 123] but in- 
cluding the possibility that amplitude of the waveform may 
also, in general, be different in an alternative theory. 

In the present work, our aim is to set up a general 
parametrization of the gravitational-wave signal in a subclass 
of ppE theories that will enable tests of GR in the radiative 
regime from GW observations, similar to the PPN and PPK 
formalisms mentioned earlier. This is an extension of our 
previous work luSll using more complete inspiral waveforms, 
called full waveforms (see Sec. lICI for a detailed discussion). 
From the ppE perspective, the model we have presented in this 
paper would correspond to the case where there are no ampli- 
tude deviations and only one phasing coefficient (correspond- 
ing to the test parameter) different from GR. In the future the 
more general ppE class of meta-models must be investigated 
to provide more generic results than those obtained from the 
subclass we deal with in this paper. 



B. Choice of PN parametrization 

The suggestion to use a smaller set of parameters comprised 
of the PN coefficient to be tested together with any two of the 
remaining PN coefficients selected as basic variables (to reex- 
press and parametrize the rest of the PN coefficients) ll25ll im- 
mediately raises the following question. Which two parame- 
ters should be chosen as the basic variables? The most natural 
choice is the two lowest-order OPN and IPN coefficients since 
they are measured most accurately. Furthermore, within GR 
at higher PN order there are spin-orbit and spin-spin terms etc, 
so that not only the choice of higher-order coefficients as basic 
parameters appears less convenient but the use of the lowest- 
order PN coefficients as basic variables may be expected to 
reduce systematic effects due to spins. Once the basic vari- 
ables are decided, the test parameter can be any of the higher 
PN order coefficients, chosen one at a time. 

The weakness of this version of the test vis-a-vis the ver- 
sion where all PN parameters are treated as independent pa- 
rameters may be worrying at first. Unlike in the latter version 
where the discrepant PN order will be explicit in the test, in 
the former version the failure of GR at a particular PN order 



may not necessarily imply the corresponding PN term to be 
different from GR. However, we believe in the robustness of 
the test itself due to the following chain of plausibility argu- 
ments: It is not unreasonable to assume that if an alternative 
theory of gravitation is consistent with GR at some PN order 
it would be normally consistent with GR at lower PN orders 
but it may disagree with GR at some higher PN order In view 
of this argument, if one is testing a particular order PN coef- 
ficient then parametrizing the lower-order PN coefficients by 
the basic PN coefficients (at OPN and IPN orders) is reason- 
able. 



What then do we hope to achieve by our choice of express- 
ing the PN coefficients of order higher than the tested PN co- 
efficient by the basic PN coefficients? To answer this ques- 
tion, we make the reasonable assumption that if a theory dif- 
fers from GR at some PN order, it is likely to differ from GR 
at higher PN orders too. Thus by parametrizing PN coeffi- 
cients of orders higher than the tested PN coefficient by the 
basic PN coefficients, we naturally take into account effects in 
GR coming from higher-order PN terms and reduce the corre- 
sponding systematic errors arising from the higher-order PN 
coefficients in the estimation of the test parameter In conse- 
quence, only departures of the correct theory of gravity from 
GR would remain and contribute dominantly to the estimate 
of the PN parameter tested in the analysis. To investigate this 
question explicitly, we computed the error in the estimation of 
a particular PN coefficient both at the 3.5PN accuracy and by 
truncating the expansion at the PN order of the test parameter. 
As expected, the two choices gave different results, with the 
full phasing yielding a more accurate estimation of parameters 
(for more details, see Sec. IIIIB3I I. Given the possibly small 
differences we are trying to explore any reduction in system- 
atic errors is indeed to be taken advantage of. 



It is interesting to note that, as pointed out in Ref. 113311 . 
the errors in the various PN coefficients we quote here can be 
translated into measurement of three- and four-graviton ver- 
tices. Keeping these caveats in view, let us consider some hy- 
pothetical theory of gravity which shows deviation from GR 
starting from 2PN order'. In our proposed test, this deviation 
would not show up when if/3 is used as test parameter, as in 
this case the deviations are only from the fact that the func- 
tional dependences of higher-order phasing coefficients on iffQ 
and i/'2 are not the same as in GR. This seems less important 
than a lower PN order test parameter itself deviating from its 
GR value. On the other hand, when iff 4 or higher PN order 
phasing coefficients are used as test parameters, this deviation 
should be evident. Thus, proceeding systematically to higher 
PN orders, one can ascertain the PN order where the new the- 
ory begins to deviate from GR. 



' There are theories which show deviation from GR starting from 2PN order. 
See for example Ref. |34il which considers one such example, though for 
spinning BHs. 



C. Scope of the current work and a summary of results 

In this work we revisit the problem in the context of the 
second generation ground-based GW interferometer such as 
Advanced LIGO and a third generation ground-based GW in- 
terferometer called Einstein Telescope (ET) that is currently 
under design study in Europe. Since ET is envisaged to have 
far better low-frequency sensitivity than Advanced LIGO (a 
lower frequency cutoff of about 1 -5 Hz), one of the aims of the 
present investigation is to evaluate the possible gains in going 
from a lower cutoff of 10 to 1 Hz. A further new ingredient in 
the present version of the above test is that we use not just the 
3.5PN RWF but also the amplitude-corrected full waveforms 
(FWF) which are 3PN accurate in amplitude (thus having 
seven harmonics other than just the leading quadrupolar one) 
and 3.5PN accurate in phase. For nonspinning binaries, the 
amplitude corrections are functions of the two masses and the 
inclination angle of the binary. The amplitude corrections at 
every PN order bring new dependences on the binary masses 
and hence could improve the estimation of the phasing coeffi- 
cients. With our previous insight in the weakening of the test 
due to the use of more parameters |'23], in the first instance, 
we make the reasonable assumption that since matched filter- 
ing procedure is more sensitive to the phase rather than ampli- 
tude, one can skip parametrizing the amplitude independently 
in the present work. Thus the deviations in the amplitude cor- 
rections of the waveform are not independently parametrized 
even though they have been taken account of in the present 
work to reduce systematic effects. We rewrite the mass de- 
pendences in the PN amplitude terms in terms of (/ro and iff2, 
just as in the case of the phase terms. 

An obvious limitation in regard to the present analysis con- 
cerns the generality of the parametrization that we employ. 
Indeed, in the strong-field regime, alternative theories may so 
qualitatively deviate from GR that the structure of the wave- 
form used here for the parametrization may not be generic 
enough to capture those features and one may need to redo 
the present analysis within a more general class of models 
like the ppE framework. It is worth mentioning that at present 
very accurate high PN order GW phasing results are available 
only for GR. Most results for alternative theories of gravita- 
tion are available only for leading or next to leading PN orders 
and hence more work is required before they are comparable 
in performance to the GR waveforms used presently for GW 
data analysis. Given these circumstances, our purpose is to 
explore what best we can extract from a subclass of PN mod- 
els close to GR and leave to future works a more complete 
investigation using e.g. ppE models. 

Based on the above analysis, within the subclass of ppE 
theories that we consider, we find that GW observations by 
advanced LIGO of binary black holes (BBHs) in the range 11- 
1 IOM0 and at a luminosity distance of 300 Mpc would allow 
the measurement of the PN coefficient i/zj, with fractional accu- 
racies better that 6% when the FWF is used. On the other hand 
using the FWF as a waveform model and a low-frequency cut- 
off of 1 Hz, observations of stellar mass BBHs in ET would 
allow the measurement of all PN parameters, except 1I/4, with 
accuracies better than 2% when the total mass of the binary is 



in the range 11-44-Mq. ET observations of intermediate mass 
BBHs would allow only two of the seven PN coefficients (1^3 
and tffsi) to be measured with fractional accuracies better than 
10% in the mass range 55-4OOM0. 

The choice of a low-frequency cutoff of 1 Hz, as compared 
to 10 Hz, reduces the relative errors in various parameters 
roughly by factors of order 2 to 10 for stellar mass black hole 
binaries. For intermediate mass binaries, which coalesce at 
smaller frequencies, though a lower cutoff helps improve the 
parameter estimation, the errors associated with the measure- 
ment of various parameters is so large that the test is not very 
interesting. Although the use of FWF has no particular ad- 
vantage in the case of stellar mass black hole binaries, their 
use in the case of intermediate mass black hole binaries im- 
proves parameter estimation by a factor of a few to almost 80. 
Large improvements are obtained for binaries that are more 
massive than ~ lOOM©. Error in the estimation of the various 
PN parameters for such systems is already so great that the 
improvement brought about by the use of FWF is not useful 
for the type of tests discussed in this study. For reasonable de- 
tection rates, intermediate mass BH binaries are at distances 
greater than 3 Gpc as opposed to 300 Mpc in the stellar mass 
case (see, e.g., Ref. 1351]): The gain due to the use of FWF is 
offset by an order-of-magnitude loss since sources are farther 
away. 

The rest of this paper is organized in the following way: In 
Sec. |II]we have introduced the noise curves that are employed 
for Advanced LIGO and ET and the waveform model used in 
the present work. This is followed by a brief description of 
the Fisher matrix formalism that will be used to perform pa- 
rameter estimation, our proposal for testing GR, the physical 
systems investigated and the implementation of the test. In 
Sec. Uni we discuss our results and various issues related to 
the present work like the systematics due to higher-order PN 
terms, choice of parametrization used in the test and the effect 
of choice of source location on the results. Finally, in Sec. |IV] 
we give a summary of our findings and future directions. 



LIGO (with SNR ~30) and observations of stellar mass (with 
SNR of few hundreds) as well as intermediate mass (with SNR 
~ 40) BBHs in ET, will allow us to measure some of the PN 
coefficients (if not all) with good accuracies. 

The first generation of long baseline interferometric 
gravitational-wave detectors (GEO600, LIGO and Virgo) 
have more or less reached their design sensitivity and have op- 
erated for a number of years taking good science quality data. 
They have shown that it is possible to build, control, and oper- 
ate highly sensitive instruments. All of these projects are now 
on the path toward building advanced versions with strain sen- 
sitivities a factor of 10 better than their current versions. This 
is made possible with research and technology in high power 
lasers, ultrahigh seismic isolation systems, improved control 
systems, etc., that has been developed over the past decade. 
When completed around 2015-2017, advanced detectors are 
expected to make routine observation of gravitational waves 
— the most promising of all sources being the coalescence 
of binaries consisting of compact objects (see Sec. Ill El for 
expected binary coalescence rates). 

While advanced detectors will open the gravitational win- 
dow for astronomical observations, the expected signal-to- 
noise ratios will not be routinely large enough to carry out 
strong-field tests of GR or high precision measurements of 
cosmological measurements. The worldwide gravitational- 
wave community has already begun to explore the technologi- 
cal development that is necessary to build detectors that are an 
order-of-magnitude better than the advanced instruments. The 
Einstein Telescope is a three-year conceptual design study 
funded by the European Commission with the goal to iden- 
tify the challenges to mitigate gravity gradient and seismic 
noise in the low-frequency region to make it possible to ob- 
serve in the 1-10 Hz band. ET will be designed to also make 
an order-of-magnitude improvement in strain sensitivity in the 
10-1000 Hz band. Such a detector will be capable of making 
routine observation of high-SNR events that will be useful for 
carrying out precisions tests of general relativity. 



II. TEST OF GR WITH ADVANCED LIGO AND THE 
EINSTEIN TELESCOPE 

In this work we shall focus on the measurement of various 
PN coefficients in the context of Advanced LIGO and the Ein- 
stein Telescope. One would like to investigate whether the 
observations of stellar mass black holes binaries in Advanced 



A. Advanced LIGO 

For the studies related to Advanced LIGO we use the Ad- 
vanced LIGO sensitivity curve IBql . The analytical fit of the 
noise curve for Advanced LIGO is given by the expression. 
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where x = f/fo, /o = 215 Hz, ^o = 10 '*'' Hz ' and /, is low- 
frequency cutoff below which 5'/,(/) can be considered infinite 



for all practical purposes. We have chosen it to be 20 Hz. The 
amplitude spectrum of Advanced LIGO is plotted in Fig.[T] 



10"' 



10"' 



^ ^10 
C/0 



10" 





'■rTT':]':;:''''!':'': ;":■;:::';!": '--Y^^J^ ■;■; ;::■;::; 


\ 


Advanced LIGO 








\ 






\, 






\ 


\ 








\ \ 


y 






N. ^^ 


^ ^ 



10' 
f(Hz) 



FIG. 1: Amplitude spectrum of Advanced LIGO and ET. 



B. Einstein Telescope 



C. The waveform model 

The earlier papers which discussed the tests of GR, includ- 
ing our own papers ll23ll25ll . assumed the so-called restricted 
post-Newtonian waveform (RWF) for quasicircular, adiabatic 
inspiral, which contains the dominant harmonic at twice the 
orbital frequency and no corrections to the amplitude. In the 
present study, we include the effects of subdominant harmon- 
ics in the waveforms. Such a wave-form, as mentioned ear- 
lier, is called the FWF and includes harmonics other than the 
dominant one, each having PN corrections to their amplitudes. 
At present the most accurate waveforms available include PN 
corrections in amplitude to 3PN order and in phase to 3.5PN 
order ll38l - l4C)ll . To see how one might test GR or, more pre- 
cisely, the structure of the PN theory, let us begin by consider- 
ing the waveform from a binary in the frequency domain. The 
full signal in its general form reads as 



The ET design study has come up with a number of possible 
sensitivity curves OTIl . Examples include a single detector 
that operates with an improved sensitivity over the whole band 
of 1 Hz to 1 kHz and a xylophone configuration consisting of 
a pair of detectors, one tuned for best low-frequency (i.e., 1- 
100 Hz) sensitivity and a second detector tuned for optimal 
performance at higher frequencies of 100 Hz to a few kHz. 
In our studies in this paper, we will use the ET-B sensitivity 
curve [331, which is also the official sensitivity curve for ET. 

An analytical fit to the ET sensitivity curve is given by 
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where x = ///o, /o = 100 Hz, Sq = lO-^^Hz"' and /, is low- 
frequency cutoff' below which 5'/,(/) can be considered infinite 
for all practical purposes. Also one has 
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The amplitude spectrum of ET is plotted in Fig.[T] In connec- 
tion with the ET design study one of the issues to be looked 
into is the science case for going down to as low a frequency 
as 1 Hz versus a more modest choice of 10 Hz. 
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Here tc and (pc are the fiducial epoch of merger and the 
phase of the signal at that epoch, respectively. Quantities 
appearing in Eq. (l2.4b with argument t{fk) denote their val- 
ues at the time when the instantaneous orbital frequency F{t) 
sweeps past the value f jk and x{t) is the PN parameter given 
by x{t) = [2nMF(t)]-^^. A(t,„/2)(f) and (f>(k,„/2)it) are the po- 
larization amplitudes and phases of the kth harmonic at «/2th 
PN order in amplitude. The coefficients in the PN expansion 
of the Fourier phase are given by 
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The constant C - 0.577 ■ ■ ■ , appearing in the expression for 
ffe, is Euler's constant. 

We have total nine post-Newtonian parameters: seven of 
these are the coefficients of v" terms for « = 0,2,3,4,5,6,7 
and two are coefficients of v"\viv terms for « = 5,6. These 
are PN coefficients in Einstein's theory and are functions of 
just two mass parameters chosen to be the total mass M and 
symmetric mass ratio v. 

In addition to mass dependence, the amplitude corrections 
also depend on the luminosity distance of the source to the 
observer and four additional angular parameters (cos ft (p, if), 
cos i) related to the source location and orientation ||4U 14211 : 6 
and determine the sources location, tfr is the polarization an- 
gle and I is the inclination angle. Once the mass dependences 
in amplitude corrections are replaced by our fundamental pair 
of i/tq and i(f2, the whole waveform can be characterized by a 
total ten parameters 
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D. Fisher matrix and statistical errors 
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We employ the Fisher matrix approach II43L 14411 to see how 
well we can measure these parameters. Below we briefly list 
the basic equations of the Fisher matrix approach that we sub- 
sequently need. 

Let ff' denote the "true values" of the parameters and let 
6" + A6" be the best-fit parameters in the presence of some 
realization of the noise. Then for large SNR, error in the esti- 
mation of parameters Aff' obeys a Gaussian probability distri- 
bution iSlfSl of die form 
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where p'"' is a normalization constant. The quantity Tab ap- 
pearing in Eq.( |2.9l l is the Fisher information matrix and is 
given by 



Fflfc = {ha I hb) 



(2.10) 



where /?„ = dhfdO". Here, ( | ) denotes the noise weighted 
inner product. Given any two functions g and h their inner 
product is defined as 



(g\h) = 4Re 



jf„„„ 



df 



nf)h{f) 

Shif) 



(2.11) 



The integration Umit {fmm,fmax\ is determined by both the de- 
tector and by the nature of the signal. Each harmonic in h{f) 
is assumed to vanish outside a certain frequency range. The 
simplest physical choice is to set the contribution from the 
A;th harmonic to the waveform zero above the frequency kf\sQ, 



where /jso is the orbital frequency at the last stable orbit 147 
Since the amplitude-corrected waveform we are using in this 
work has eight harmonics, we set the upper cutoff to be 8/iso 
when we use the FWF in the analysis. For lower cutoff, as 
power spectral densities S i,(f) tend to rise very quickly below 
a certain frequency fs where they can be considered infinite 
for all practical purposes, we may set it to be fs- Using the 
definition of the inner product one can reexpress T„h more ex- 
plicitly as 



ab 



rkfw 



Re(h:(f)hb(f)) 
Shif) 



df. 



(2.12) 



The covariance matrix, defined as the inverse of the Fisher 
matrix, is given by 



r'' s (Afl^Afl'^) = (F-')"*, 



(2.13) 



where (■) denotes an average over the probability distribution 
function in Eq. ( 12. 9t . The root-mean-square error cr^ in the 
estimation of the parameters ff' is 



<(Aft')^> 



2\l/2 



(2.14) 



In the present work we deal with inspiralling compact bina- 
ries as seen by Earth bound detectors. For such burst sources, 
one can approximate the detector's beam pattern functions 
as being constant over the duration of the signal and thus 
we can assume that angular parameters (cos ft <p and ifr) as 
well as the luminosity distance (Di) are fixed and thus can 
be excluded from the analysis. With this restriction, the 
large 10-dimensional parameter space reduces to a smaller 6- 
dimensional parameter space given by 



p = (fc, 4>c, lAo, lAa, lAr, cos i) 



(2.15) 



In order to test the PN structure, one should be able to 
measure various PN coefficients with good accuracy. In the 
present work we have assumed that the relative error in the 
measurement of a parameter should be less than 10%, i.e. 
Aif/jfiffj < 0.1, where Aif/j is the error in the estimation of the 
parameter if/j, in order to estimate its value in the PN series 
with confidence. 



E. Systems investigated 

The first detection of gravitational radiation in ground- 
based interferometric detectors is generally expected to be 
from the coalescence of compact binary systems with neutron 
star and black hole components ll48ll . Among these, binary 
neutron stars (BNS) are arguably the most promising ones 
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FIG. 2: Plots showing the regions in the irii-mi plane that correspond 
to l-o" uncertainties in ifio, i//2 and tf/s] (left panel) and those in if/Q, 
tf/2 and i/'5imod (right panel) for a (2, 20)Mq BBH at a luminosity dis- 
tance of Di = 300 Mpc observed by ET. The low-frequency cutoff 
is IHz and RWF has been used. For the curves in the right panel 
we have assumed that the correct theory of gravity is a hypothetical 
non-GR theory in which the phasing coefficient tf/si and all higher PN 
coefficients differ from the GR values by 1%. 



with expected rates of about 40 mergers per year in Advanced 
LIGO and millions of them in ET. While very interesting for 
other proposed tests of GR, BNS systems are not useful for 
the tests proposed in this study. For our purposes a compact 
binary in which one or both the components is a stellar mass 
(~ 2-3OM0) or intermediate mass (~ SO-IOOOMq) black hole 
(the other being a neutron star) would be most interesting. For 
our studies related to Advanced LIGO, we have chosen binary 
black holes in the mass range II-IIOM0 and their distance 
from the Earth to be 300 Mpc. 



For the analysis using ET we have discussed separately 
stellar-mass and intermediate-mass BBHs. For stellar mass 
BBHs, we have again chosen their luminosity distance from 
the Earth to be 300 Mpc and the range of the total mass to be 
1 l-44Mo. Coalescence rate of stellar mass BBHs is highly un- 
certain. The predicted rate of coalescence within a distance of 
300 Mpc varies between one event per 10 years to several per 
year |J49[]. However, it is with such rare high-SNR events that 
one expects to perform precision tests of GR. For intermedi- 
ate mass black holes, we have chosen the distance to be 3 Gpc 
(z = 0.55), and their total mass to be in the range 55-11OOM0. 
The evolutionary history of intermediate mass BBHs and their 
rate of coalescence is still not well understood. The main mo- 
tivation to study these systems comes from the models that in- 
voke them as seeds of massive black holes at galactic nuclei. 
In a recent study, it has been suggested that only few coales- 
cence events of intermediate mass BBHs could be expected 
within a redshift of z -2. Also depending on what triggered 
seed galaxies there may be a few events within a redshift of 

z = 1 illisHIl. 



F. Implementation of the test 



As mentioned earlier, in Einstein's theory (and thus in the- 
ories "close" to GR) each PN coefficient for a nonspinning 
compact binary is a function of the two mass parameters, the 
total mass M and the symmetric mass ratio v. In other words, 
we can say that each (/r, is a function of the masses (oti, 1112) of 
the components constituting the binary, i.e. i/r, s i^,(mi,m2). 
With high-SNR GW observations of stellar and intermediate 
mass BBHs in Advanced LIGO and ET, it would be possible 
to measure the individual masses constituting the binary with 
good accuracies. Thus, once the (statistical) error in the pa- 
rameter is estimated using the Fisher matrix, we can represent 
the region it spans in the space of masses by inverting the re- 
lation 1^/ s t/^i(mi,m2) to get say m2 = m2(i//i,mi). Given the 
measured value i/r™*^"^ and the errors Ai/f,- in the estimation of 
(/',, the region in the mass-plane corresponding to 1112 is given 
by 1112 = m2(i^™'^'"* ± Aiffi, mi). For each ^i, there would be an 
allowed region in the 1111-1112 plane and if Einstein's theory of 
gravity, or, more precisely, the PN approximation to it, is a 
coiTect theory then the three parameters 1^0, 1A2 and 1^7- (the 
test parameter) should have a common nonempty intersection 
in the mi-m2 plane. Proceeding in this way, for six test pa- 
rameters we shall have six different tests of the theory. In the 
present work, we shall only discuss asymmetric binaries with 
component mass ratio q,„ =0.1. Since the different PN coef- 
ficients are symmetric with respect to the exchange of mi and 
1112, we expect plots in the 1111-1112 plane to have two symmetric 
branches. Figures |2]|6] and [8]show one branch of the full plot. 

Figure |2] schematically demonstrates how the test works by 
using iff() and 1^2 as basic variables and ^51 as a test param- 
eter. The plot on the left uses PN coefficients predicted by 
GR, assuming GR is a correct theory of gravity. Clearly, this 
shows that all three parameters, iffQ, ^2 and ^si, have a com- 
mon nonempty intersection in the plane of masses and this is 
what we expect if GR is the correct theory of gravity. 

In contrast, consider the possibility that the correct theory 
of gravity is a hypothetical non-GR theory in which the phas- 
ing coefficient 1^5/ and all higher PN coefficients, differ from 
the GR values by 1%. We have assumed here, in an ad-hoc 
manner, that the deviation of the PN terms at higher orders 
above the 2.5PN term (which is put to test) to be a simple 
scaling, i.e., i/r^t— >l-01i/'jt for A: > 5. As we shall show later, 
there is a range of binary masses for which the observation of 
the GW signal by ET could estimate this coefficient with an 
accuracy much better than 1%. In this scenario, if we inter- 
pret the i/'5i obtained by fitting to the observed GW signal, as 
a GR coefficient, there will definitely be an inconsistency in 
the mi -m2 plane. This can clearly be seen in the right panel of 
Fig. |2] where there is no overlapping region between the three 
parameters in question in the 1111-1112 plane, thus demonstrating 
the spirit of the proposed test. 

Figure[3]shows a similar exercise for the (10, 1OO)M0 BBH 
located at a luminosity distance of Dl = 3Gpc observed by 
ET. The low-frequency cutoff is IHz and RWF has been used. 
The left panel in Fig. [3] assumes that GR is the correct theory 
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FIG. 3: Plots showing the regions in the nii-mi plane that correspond to 1-cr uncertainties in Newtonian, IPN and 2.5PN coefficients in the PN 
series for a (10, 100)Mo BBH at a luminosity distance of D/, = 3 Gpc observed by ET. The low-frequency cutoff is IHz and RWF has been 
used. The left panel corresponds to OR as the correct theory of gravity while the middle and right panels correspond to hypothetical non-GR 
theories of gravity which have phasing coefficients xj/si (2.5PN) and higher differing from the GR values by 1% and 10% respectively. 



of gravity, whereas the middle and right panels assume that 
the correct theory of gravity is a hypothetical non-GR theory 
in which the PN coefficients at 2.5PN (i.e., 1^5;) and all higher 
orders differ from their GR values by 1% and 10% respec- 
tively. The 2.5PN coefficient in GR and in the above two hy- 
pothetical theories can be measured with fractional accuracies 
of 5.6%, 5.5% and 5.1%, respectively, for the system under 
consideration. 

As in the previous case (cf.. Fig |2|, we notice the depar- 
tures from GR. However, these departures are not very clear 
in the middle panel where the 2.5PN coefficient and other 
higher-order terms differ from their GR values by only 1%. 
As a result one would not be able to discriminate between the 
two theories by using GW observations of such sources in ET, 
even though the parameter can be measured accurately. The 
right panel of the Fig. [3j which corresponds to a theory in 
which the values of 2.5PN and higher-order coefficients differ 
from their GR values by 10%, brings such significant depar- 
tures from GR that despite the larger errors in the estimation 
of the test parameter, one can distinguish between the two the- 
ories using the test we are proposing here. One should bear in 
mind that the model we have used is the simplest possible way 
in which a deviation could occur. But our purpose here is to 
have a proof of principle demonstration of the power of the 
proposed test, given the lack of knowledge of the exact man- 
ner in which such deviations could occur 

To summarize, assuming that GR is correct, our proposal 
gives the accuracy with which three of the PN parameters can 
be measured. How does that test GR? For e.g. if GR is not 
correct and differs at, say 1 .5 PN level onwards, then our claim 
is that we would begin to see inconsistencies in the estimated 
parameter values beyond the accuracy of measurement pro- 
vided deviations from GR are large enough. One may be con- 



cerned about the extent to which the departure of higher-order 
terms from their GR values would penalize the estimation of 
lower-order terms. As evidenced by our examples above, they 
induce bias in the estimation of parameters but do not lead to 
greater errors in the estimation of parameters. In other words, 
the 1.5PN and higher-order PN coefficients not agreeing with 
GR might shift the mean of the distribution of (M, v) but the 
width should remain more or less the same. Put differently, 
if the PN expansion differs from GR slightly then the error in 
the estimation of parameters will not change to first order. 



III. THE RESULTS 
A. Advanced LIGO 

In this section we investigate the possibility of the test using 
GW observations of BBHs in Advanced LIGO. As discussed 
earlier, the range of total mass explored is 11-110 Mq and we 
assume that binaries are located at a luminosity distance of 
300 Mpc. Plots in Fig.|4]show the variation of relative accura- 
cies with which two of the PN coefficients, 1^3 and i/'si, can be 
measured by Advanced LIGO using the restricted and the full 
waveform. The components of the binary have the mass ratio 
of 0.1. 

It is evident from the plots that when the FWF is used, 1/^3 
and 1^51 can be measured with fractional accuracies better than 
6% and 23%, respectively, in the whole mass range under con- 
sideration. We shall require (rather arbitrarily) that the relative 
error in the measurement of a PN coefficient be less than 10% 
in order for the test to be effective. Clearly, 1^3 can be es- 
timated quite accurately and thus it can be used to test the 
theory. On the other hand, since 1^51 is not so well determined. 
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FIG. 4: Plots showing the variation of relative errors t^tfjjjxfjj in the test parameters i/'j- = i/'s, 1^5/ as a function of total mass of binaries in the 
range 11-110Mq (with component masses having mass ratio of 0.1) located at 300 Mpc observed by Advanced LIGO, using both the RWF 
and the FWF as waveform model with the source orientations chosen arbitrarily to be = = 7r/6, ii = n/4, and i = n/3. The noise curve 
corresponds to the one shown in Fig.[T]for the Advanced LIGO case and its analytical fit is given by Ea.l l2.lt . It is evident from the plot in the 
left panel that the fractional accuracies with which 1(13 can be measured are better than 6% for the entire mass range under consideration when 
FWF is used and thus can be used to test the theory of gravity, tf/st (right panel) can be measured with fractional accuracies better than 23% 
for the entire mass range when FWF is used but being a poorly determined parameter it can provide a much less stringent test of the theory of 
gravity. 



it can still provide a less stringent test of the theory. The mea- 
surement of other PN coefficients is not accurate enough to 
lead to a meaningful test of GR. 

The plots clearly show the benefits of bringing higher har- 
monics into the analysis. The use of the FWF typically im- 
proves the estimation by a factor of 3 to almost 100. 



B. Einstein Telescope 

In the previous section we have seen that with Advanced 
LIGO one can only test PN theory up to 1.5PN. Can one do 
better with the proposed third generation detector like the ET? 
In what follows we investigate the extent to which one can test 
the PN theory using GW observations of stellar mass and in- 
termediate mass BBHs using ET. In addition to this we will 
discuss some other key issues influencing the results such as 
eff'ects of PN systematics on the test, choice of parametriza- 
tion and dependence of the test on angular parameters. 



10 vs 1 Hz. The top and bottom panels correspond to the 
lower frequency cutoff of 1 and 10 Hz, respectively, while 
the left and right panels correspond to the RWF and FWF, re- 
spectively. The source orientations are chosen arbitrarily to 
he - (p - n/6, 41 - n/A, and l - n/3. It should be evi- 
dent from the plots that the best estimates of various test pa- 
rameters are for the combination using the FWF with a lower 
cutoff frequency of 1 Hz. In this case, all i^,'s except 1^4 can 
be measured with fractional accuracies better that 2% for the 
total mass in the range 11-4-4-Mq. On the other hand when the 
lower cutoff is 10 Hz, with the FWF all i^,'s except 1^4 can be 
measured with fractional accuracies better than 7%. It is also 
evident from the plots that as compared to other test parame- 
ters, 1^3 is the most accurately measured parameter in all cases 
and best estimated when the lower frequency cutoff is 1 Hz. 
On the other hand, 1/^4 is the worst measured parameter of all 
the test parameters. However, we see the best improvement in 
its measurement when going from the RWF to the FWF. 



1. Stellar mass black-hole binaries 

Figure |5] plots the relative errors Aij/tIiPt as a function of 
total mass M of the binary at a distance of Di - 300 Mpc. 
We have considered stellar mass BBHs of unequal masses and 
mass ratio 0.1, with the total mass in the range ll-44Mo. Fig- 
ure|5]also shows two types of comparisons: (a) full waveform 
vs restricted waveform, and (b) a lower frequency cutoff of 



Figure |6] shows the regions in the nii-nio plane that cor- 
respond to 1-cr uncertainties in (//q, ^2 and various test pa- 
rameters which in turn will be one of the six test parameters 
4iT - iA3>'A4''A5;>'A6-iA6/, and 1/^7, one at a time, for a (2, 20) 
Mq BBH, at a luminosity distance of Dl - 300 Mpc observed 
by ET. It is evident from the plots coiTesponding to various 
tests that each test parameter is consistent with coiTesponding 
fundamental pair (tf/Q, if/o). 
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FIG. 5: Plots showing the variation of relative errors Atf/T/if/T in the test parameters iffT=tf/3, tf/4, if/51, ifre, 4'e>h and 1/^7 as a function of total mass 
M for stellar mass black hole binaries (with component masses having mass ratio 0.1) at a luminosity distance of D^ = 300 Mpc observed by 
ET, using both RWF (left panels) and FWF (right panels) as waveform models. The choice of the source orientations is the same as quoted 
in Fig. |4] The noise curve corresponds to the recent ET-B sensitivity curve. Top panels correspond to the lower frequency cutoff of 1 Hz. 
By using FWF as the waveform model all i/'^t's except ^4^4 can be tested with fractional accuracy better than 2% in the mass range ll-44Mo. 
Bottom panels correspond to the lower frequency cutoff of 10 Hz. Using FWF, all i^^'s except 1/^4 can be tested with fractional accuracy better 
than 7% in the mass range 1 l-44Mo. 



Intermediate mass black hole binaries 



and L - nITi. 



Figure I2] plots the relative errors Aij/j Ixjjt as a function of 
the total mass M of the binary at a distance of Di-1> Gpc. We 
have considered BBH of unequal masses with mass ratio 0. 1 . 
As in Fig.|5] Fig.|2]also shows two types of comparisons: (a) 
the effect of the use of FWF on parameter estimation against 
RWF, (b) the effect of lowering the cutoff frequency from 10 
to 1 Hz. As before, top and bottom panels correspond to the 
cutoff frequency of 1 and 10 Hz, respectively, and left and 
right panels to RWF and FWF, respectively. The source ori- 
entations are chosen arbitrarily ioht 9 - (p - n/6, i// = 7r/4, 



It is evident from the plots that the least relative errors in 
various test parameters are for the combination that uses the 
FWF and a lower cutoff of 1 Hz. Unlike the case of stel- 
lar mass BBHs, in the case of intermediate mass BBHs only 
two of the test parameters, tfr^ and (/f^/, can be measured with 
fractional accuracies better that 10% for the total mass in the 
range 55-400 Mq with FWF and lower cutoff frequency as 1 
Hz. On the other hand, when the lower frequency cutoff is 
10 Hz the use of the FWF allows the estimation of 1/^3 and 1^5; 
with fractional accuracies better than 10% for the total mass in 
the range 90-220 Mq. As compared to other test parameters. 
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FIG. 6: Plots showing the regions in the 1121-1112 plane that correspond to l-cr uncertainties in tf/o, 1^2 and various test parameters, which happen 
to be one of the six test parameters i/tj = ifi3,^4,il/^i,i//^,^(,i,andi//j atone time, for a (2, 20) Mq BBH at a luminosity distance of £)/. = 300 Mpc 
observed by ET. In all six plots shown above t//o and i/f2 are chosen as the fundamental parameters (from which we can measure the masses of 
the two black holes). Each parameter corresponds to a given region in the nii-/)j2 -plane and if GR is the correct theory of gravity then all three 
parameters, tf/Q, 1/^2 and tf/j should have a nonempty intersection in the mi-nu plane. A smaller region leads to a stronger test. Notice that all 
panels have the same scaling except the top middle panel in which Y axis has been scaled by a factor 10. 



1^3 is the most accurately measured parameter in all cases and 
best estimated when the low-frequency cutoff is 1 Hz. Param- 
eters i//4 and i//(, are poorly measured as compared to the other 
test parameters but again we see the best improvement in the 
estimate of t/r4 when using the FWF. 

Figure [8] shows the regions in the m\-m2 plane that coiTe- 
spond to l-cr uncertainties in if/Q, ^2 and the test parameters 
if/T - iA3' ^A, ^1^51, lAe- ^bi, tAv^ one at a time, for a (20, 200) Mq 
BBH at a luminosity distance of Di-1> Gpc observed by ET. 
It is clear from the plots that each test parameter is consistent 
with the corresponding fundamental pair (i/tq, i^a)- 



3. Effects ofPN systematics on the test 

The inability to measure all the PN parameters simultane- 
ously led us to propose a more modest procedure to test the 
PN parameters one at a time. In parameter estimation, it seems 
intuitive not to ignore our knowledge of the known high PN 
order phasing. Further, it is natural to assume that if an al- 
ternative theory of gravitation, similar to GR, agrees with GR 
at some PN order, it would agree with it at a lower PN order 



but may differ from it at some higher PN order. Thus, when 
testing a coefficient at some particular PN order, expressing 
a lower-order PN coefficient in terms of the basic pair of PN 
variables seems reasonable. However, expressing the higher 
PN order coefficients in terms of the basic pair may appear 
more disconcerting. Here we look at the issue in a little more 
detail and provide our point of view. 

We propose a comparison of the following two schemes: 
The first scheme, as before, uses if/o and )/'2 as basic parameters 
and all known PN coefficients up to 3.5PN, except the test 
parameter, are expressed in terms of if/o and )/'2. The second 
scheme is similar but the phase evolution is truncated at the 
PN order corresponding to the test parameter. 

Thus, to test i/'3 in the second scheme the phasing is trun- 
cated at 1.5PN, to test i/'4 at 2PN and so on. Figure |9] com- 
pares the two schemes. It should be evident from the figure 
that in the first scheme the use of 3.5PN phasing, rather than 
a lower PN order ( e.g., 1.5PN in testing xj/^), does improve 
the accuracy with which one can measure a certain parame- 
ter. Conversely, the poorer estimate ( i.e. larger error) in the 
second scheme is due to the neglect of higher PN order terms. 

One can infer, therefore, that what is achieved in the first 
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FIG. 7: Same as Fig|5]but for intermediate mass black hole binaries (witli component masses liaving mass ratio 0. 1) at a luminosity distance of 
£)l=3 Gpc. With lower frequency cutoff of 1 Hz, using FWF as the waveform model, i/'3 and i/'s/ can be tested with fractional accuracy better 
than 10% for the mass range 55-400Mo. On the other hand, with a lower frequency cutoff of 10 Hz, using the FWF, \fj^ and i/^s/ can be tested 
with fractional accuracy better than 10% for the mass range 90-220Mo. 



scheme is an improvement in parameter estimation arising 
from higher-order PN phasing; parametrization of higher- 
order PN coefficients in terms the two basic variables is in- 
deed a reasonable choice if one wants to quantitatively look 
at the deviation from standard GR of an alternative theory of 
gravitation similar to GR^. Let us also note that if a theory 
of gravity deviates from GR at a particular PN order, then in 
our first scheme the test may actually fail at a lower PN or- 



der As a result, our test will not be able to conclusively assert 
the PN order that is inconsistent with GR. Rather a failure of 
our test is indicative of the failure of GR at some PN order. It 
would then be necessary to carry out a more powerful test of 
the theory by treating all the PN coefficients as independent 
parameters. Such a test could in principle help determine at 
which order(s) the true theory of gravity is inconsistent with 
GR. 



^ A variant of the test starting from GW phasing expressions in a more gen- 
eral meta theory may need to be implemented to include theories with dif- 
ferent PN structure like scalar-tensor theories with qualitatively different 
effects like dipole radiation. 



4. The choice of basic parametrization and the accuracy of the test 

As discussed in Sec.|l] in the present work we have chosen 
the lowest-order (and hence the best determined) PN coeffi- 
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FIG. 8: Same as in Fig. |6]but for intermediate black hole binaries in the mass range (20, 200)Mo at a luminosity distance of J9/,=3Gpc 
observed by FT. Notice that all bottom panels and the top middle panel have the same scaling whereas the Y axes of the top left panel and the 
top right panel have been scaled by a factor of 5. Note that there appears just one boundary for 1/^4 in the plot shown in top middle panel since 
the other bound does not exist for the range of values on X axis. 



cients (/tq and 4/2 to parametrize the waveform. One might 
wonder whether the choice of i/'o and ^2 as basic variables 
is the most optimal. To investigate this further, we explored 
other choices of the basic pair to parametrize the waveform, 
e.g. {if)(),if/f) where / can be one of (3, 4,5/, 6, 61, 7). Ta- 
ble U shows a comparison of the accuracies of the measure- 
ment of the various PN parameters under different choices of 
the parametrization schemes for a (10, 1OO)M0 binary, located 
at a luminosity distance of 3 Gpc. The comparison uses the 
RWF, with lower frequency cutoffs of 1 and 10 Hz, respec- 
tively. 

From the table the following observations are evident: 

1 . A comparison of the values in blocks symmetric across 
the principal diagonal one can compare the errors in the 
estimation of a particular parameter in the following 
two cases: once when the parameter is one of the ba- 
sic variables and secondly when it is a test variable. It 
is also clear that, in general, a parameter is determined 
more precisely when it is a basic parameter than when 
it is a test parameter. This is mainly because the basic 
variables bring new functional dependences via the rest 
of the phasing terms. 



2. The choice of the lower-order PN coefficient (/'2 as a 
basic variable leads to a more precise test. 

3. When iItq is one of the basic variables, the dispersion 
in the relative error of the other basic variable is least 
when the lowest-order PN coefficient t/'2 is chosen as 
the second basic variable, 

4. An interesting case corresponds to the choice of 1/^4 as 
the basic variable which seems to allow for the best de- 
termination of i//4. 

As a result, although, in principle, one has the freedom of 
parametrizing the waveform in terms of any of the two PN 
coefficients, the choice of (i/^o , 1A2) as basic variables is the 
optimal one. 

The above question may be equivalently investigated by 
looking at the volume of the three-dimensional ellipsoid cor- 
responding to the three phasing coefficients which are in- 
volved in the test. We find that the smallest volume of the 
ellipsoid corresponds to the case where i/^o and if)2 are used 
as basic variables as compared to other combinations, for all 
test parameters except 1/^4 and 1/^6. For these two parameters 
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FIG. 9: A comparison of relative errors in the measurement of various PN parameters for binaries with masses in the range 55-220Mo at a 
luminosity distance of 3 Gpc for two cases: The first is the same as before when i/'o and \fj2 are basic parameters and all PN parameters up to 
3.5PN (full phasing) except the test parameter are parametrized by i/'o and i/'2- The others are similarly constructed but the phasing is truncated 
at the PN order corresponding to the test parameter. The low-frequency cutoff is 10 Hz and the RWF has been used. The test parameter with 
truncated phasing is denoted by i/r„ while with full 3.5PN phasing it is denoted by t/^/f . 



the volume is smaller when they are used as basic variables 
together with i/'o- 



5. The choice of angles 

In SeclUwe pointed out that the signal depends on four an- 
gular parameters (cos Q, 0, i^, cos i) related to the source loca- 
tion and orientation but that they were chosen arbitrarily to be 
6 - (p - 7t/6, t// - njA, and l - njli in the present study. This 



is because for terrestrial detectors and burst sources the angles 
could be considered constant. To quantify the effect of these 
angular dependences on the test we computed the relative er- 
ror in a particular PN parameter for 100 different realizations 
of these angular parameters. 

The result is plotted in Fig. [TOl From Figs. [TOl and |7] it is 
clear that the value of the relative error in the estimation of i/'3 
for a (10, 1OO)M0 binary located at a luminosity distance of 
3 Gpc is a typical value and, as was physically expected, the 
weak dependence on angles is a good approximation. 
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TABLE I: Accuracies of the measurement of various PN parameters (using tlie RWF and 1 Hz low-frequency cutoff) for (10, 100)Mo binary 
located at a luminosity distance of 3 Gpc, with different choices of parametrization schemes. For each entry, the number within parentheses is 
the factor by which the accuracy will be reduced if a lower cutoff of 10 Hz is chosen instead of 1 Hz. For the fundamental pair we have chosen 
(i^o> lA/) where / can be any of 2, 3, 4, 5/, 6, 61, 7. In each case, the relative error in the test parameter is listed in the third row. 
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IV. SUMMARY AND CONCLUDING REMARKS 

In this paper we have studied the possibility of testing the 
theory of gravity within a well-defined subclass of ppE theo- 
ries using GW observations of BBHs by a typical second gen- 
eration GW interferometer (Advanced LIGO) and the plausi- 
ble third generation GW interferometer (ET). Within this sub- 
class of theories for Advanced LIGO we have shown that GW 
observations of BBHs (in the range 11-110Mq and at a lu- 
minosity distance of 300 Mpc) can be used to estimate only 
the PN coefficient i^s with fractional accuracy better than 6% 
when the FWF is used (see Fig. HJi. Estimation of a PN co- 
efficient with such an accuracy suggests that Advanced LIGO 
could indeed begin the era of strong-field tests of gravity. We 
have also compared the results for the FWF and RWF and 
shown that FWF reduces the errors by a factor of 3 to almost 
100. 

We have also studied in detail the stellar mass and inter- 
mediate mass regimes of the compact binary source popula- 
tion in the ET sensitivity band, for 1 and 10 Hz lower cut off 
frequencies and compared the advantage of using the FWF 
model over the RWF model. We find that the lower frequency 
cut off of 1 Hz plays a crucial role in testing GR with ET. 

For stellar mass binary coalescences (total mass < 44M0) 
as well as intermediate BH binaries, the lower cut off of 1 Hz 
improves the estimation of all PN parameters in the phasing 



formula. For stellar mass binaries, the improvement in the es- 
timation is between a factor of 2 to almost 20 when the RWF 
model is used. When the FWF model is used, the improve- 
ments are typically between factors 2-10 (see Fig.|5]l. 

For intermediate mass binaries, which coalesce at lower fre- 
quencies, though the smaller lower cut off improves the pa- 
rameter estimation, the errors associated with the measure- 
ment of various parameters are so large that the tests are not 
very interesting (see Fig. |7]). However, when total mass is 
less than about lOOM©, all the i/r^^'s are measured with relative 
errors less than unity, the most accurately determined param- 
eters being i/'3 and 1^5/, which are determined with accuracies 
better than 10%. This seems to be the most interesting mass 
range for the proposed test in the ET band. Though the use 
of the FWF does improve the estimation of various parame- 
ters, the test is less impressive since for astrophysically real- 
istic event rates, we have to consider distances as large as 3 
Gpc (as opposed to 300 Mpc for the stellar mass case). Thus, 
only if there is such an event very close by, can the test be 
performed very accurately. 

It is worth bearing in mind that in addition to systematic 
effects due to higher-order PN terms, various other systematic 
effects could offset the accuracy of the proposed test: 

1 . If the components of the binaries have spins, the phas- 
ing coefficients are functions not only of the individual 
masses but also the spin parameters. Further, if the bi- 
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nary is precessing (which would be the case when the 
spins are not aligned or antialigned with the orbital an- 
gular momentum vector of the binary), the waveforms 
will have a very different structure due to spin-induced 
modulations. To get a simple estimate of the effect of 
spins on the proposed test, we consider spinning but 
nonprecessing binaries. For such binaries, the effect 
of spin is to introduce additional spin-dependent con- 
tributions in various phasing coefficients at and above 
1.5PN. The 1.5PN phasing coefficient in this case has 
an additional spin parameter [}, which is a function 
of the individual spins of the binary and takes values 
< /? < 8.5 (sl. We found that for values of /? > 6, 
the bias in our estimate could be more than 100%. This 
means that the presence of spins could significantly bias 
the proposed test for large values of the spin parameter. 

Another effect is that of the orbital eccentricity, which 
we have ignored by assuming the binary's orbit to be 
quasicircular As shown in Refs. jlq, 15311 . orbital ec- 
centricity will introduce additional phasing coefficients 
with completely different frequency dependences. It 
will need a careful study to assess how to incorporate 
the effect of eccentricity into our analysis, which we 
postpone to a future work. 

Lastly, since we have used PN inspiral waveforms, the 
neglect of merger and ringdown effects could also lead 
to further systematic errors. By a proper choice of the 
domain of integration of the signal, we should be able 
to take care of it to some extent. A detailed study us- 
ing some of the analytic parametrizations of numerical 
relativity waveforms (see Refs. ll54l - t56ll ) is planned as a 
follow up of this work. 



The extension of these results within a more general class of 
models like the ppE-framework can be expected to provide 
more general results in the future. 

To fully test our proposal one must mimic the whole exer- 
cise with mock data. One has to inject a non-GR signal into 
Gaussian background with a signal that differs from GR at 1 .5 
PN and higher orders by a certain degree. One would then 
need to extract, say, the first three parameters by an Markov 
chain Monte Carlo technique that employs GR templates, and 
see if what we expect based on our toy examples above holds 
good. We are currently exploring this exercise but this exer- 
cise is quite compute intensive and goes beyond the scope of 
the present paper and so we defer its full discussion to a future 
publication. 




FIG. 1 1 : This plot shows the variation of systematic bias due to spin 
F{fi) with the spin parameter < /3 < 8.5, where F{fj) is given by 
F(je) = 4/3(16;r-4/3r'. 



60 
50 



>,40h 
o 

CD 

O" 
CD 



20- 



10 



0.02 



I c 

0.04 

AV3/V3 



0.06 



0.08 



FIG. 10: Histogram for the relative error in the estimation of the 
parameter 1^3 using 100 different realizations of angular parameters 
for a (10, 100)Mo binary located at the luminosity distance of 3 Gpc. 
The low-frequency cutoff is 1 Hz and RWF has been used. 

In this paper we have explored issues and indicated ways to 
test a class of theories of gravitation close to GR by using GW 
observations in Advanced LIGO and the Einstein Telescope. 
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Appendix A: Systematic effect due to spin 

We discuss the typical biases on our estimates due to the 
assumption that the binary components are nonspinning. We 
demonstrate this, by taking the 1.5PN phasing coefficient, 
where the spins first enter the phasing. For convenience, we 
have assumed the spins of the binary are aligned with the or- 
bital angular momentum vector, in which case we can use the 
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direct analytical formula for the phasing coefficient. 

As we mentioned earlier, the nonspinning 1.5PN phasing 
coefficient is given by q,"°"''P'" - -I5;r. The corresponding 

expression for spinning but nonprecessing binaries is aj'"" = 
-I6n + 4/3 where yS is a spin parameter which is a function 
of the spins of the individual components of the binary and 
whose value lies in the range < /? < 8.5 ll52ll . Thus the 
difference in the value of the 1.5PN coefficient due to spin is 
- 4/3. The bias in our estimates of 
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where F(J3) - ^g^ 4^-, quantifies the bias in our estimate. 

Figure [TT] shows the plot of F(J3). As is obvious, the sys- 
tematic bias due to spins could offset the estimation of ot, by 
more than 100% for/3 > 6. 



[1] C. M. Will, Living Rev. Rel. 9, 3 (2006), gr-qc/05 10072. 

[2] C. M. Will, Theory and experiments in gravitational physics 

(Cambridge University Press, New York, USA, 1981). 
[3] T. Damour and J. H. Taylor, Phys. Rev. D 45, 1840 (1992). 
[4] T Damour and N. Deruelle, Ann. Inst. H. Poincare Phys. Theor. 

43, 107 (1985). 

[5] T. Damour and N. Deruelle, Ann. Inst. H. Poincare Phys. Theor. 

44, 263 (1986). 

[6] J. Taylor and J. Weisberg, Astrophys. J. 253, 908 (1982). 

[7] B. S. Sathyaprakash and B. F. Schutz, Living Rev. Rel. 12, 2 

(2009), arXiv:0903.0338. 
[8] A. Buonanno, B. R. Iyer, E. Ochsner, Y. Pan, and B. S. 

Sathyaprakash, Phys. Rev. D 80, 084 043 (2009), 0907.0700. 
[9] |http : //www ■ ligo. caltech . edu| 
[10] 'http : //www . Virgo . infn . it ' 
[II] L. Blanchet and B. S. Sathyaprakash, Class. Quantum Grav. 11, 

2807 (1994). 
[12] L. Blanchet and B. S. Sathyaprakash, Phys. Rev. Lett. 74, 1067 

(1995). 
[13] F. D. Ryan, Phys. Rev. D 56, 1845 (1997). 
[14] C. M. Will, Phys. Rev. D 50, 6058 (1994), gr-qc/9406022. 
[15] C. M. Will, Phys. Rev. D 57, 2061 (1998), gr-qc/970901I. 
[16] A. Krolak, K. D. Kokkotas, and G. Schafer, Phys. Rev. D 52, 

2089 (1995). 
[17] C. M. Will and N. Yunes, Class. Quantum Grav. 21, 4367 

(2004), gr-qc/0403 100. 
[18] E. Berti, A. Buonanno, and C. M. Will, Phys. Rev. D 71, 084025 

(2005), gr-qc/041I129. 
[19] K. G. Arun and C. M. Will, Class. Quant. Grav. 26, 155002 

(2009), arXiv: 0904.1190. 
[20] A. Stavridis and C. M. Will, Phys. Rev. D 80, 044002 (2009), 

arXiv:0906.3602. 
[21] K. Yagi and T Tanaka, Phys. Rev. D 81, 064008 (2010), 

0906.4269. 
[22] L. S. Finn and P J. Sutton, Phys. Rev. D 65, 044022 (2002). 
[23] K. G. Arun, B. R. Iyer, M. S. S. Qusailah, and B. S. 

Sathyaprakash, Class. Quantum Grav. 23, L37 (2006), gr- 

qc/06040I8. 
[24] M. S. S. Qusailah, Phd thesis, Jawaharlal Nehru University, 

New Delhi (2006). 
[25] K. G. Arun, B. R. Iyer, M. S. S. Qusailah, and B. S. 

Sathyaprakash, Phys. Rev. D 74, 024006 (2006), gr- 

qc/0604067. 
[26] O. Dreyer, B. Kelly, B. Krishnan, L. S. Finn, D. Garrison, and 

R. Lopez-Aleman, Class. Quantum Grav. 21, 787 (2004), gr- 

qc/0309007. 
[27] E. Berti, V. Cardoso, and C. M. Will, Phys. Rev. D 73, 064030 

(2006), gr-qc/05 12160. 
[28] S. A. Hughes and K. Menou, Astrophys. J. 623, 689 (2005), 



astro-ph/04I0I48. 
[29] D. Keppel and P Ajith (2010), arXiv: 1004.0284. 
[30] S. Alexander, L. S. Finn, and N. Yunes, Phys. Rev. D 78, 

066005 (2008), 0712.2542. 
[31] C. Molina, P. Pani, V. Cai'doso, and L. Gualtieri Phys. Rev. 

D 81, 124021 (2010), arXiv:I004.4007[gr-qc]. 
[32] N. Yunes and F Pretorius, Phys. Rev. D 80, 122003 (2009), 

0909.3328. 
[33] U. Cannella, S. Foffa, M. Maggiore, H. Sanctuary, and R. Stu- 

rani, Phys. Rev. D 80, 124035 (2009), 0907.2186. 
[34] C. F Sopuerta and N. Yunes, Phys. Rev. D 80, 064006 (2009), 

arXiv:0904.4501. 
[35] I. Mandel, J. R. Gair, and M. C. Miller (2009), 

arXiv:0912.4925. 

[36] Tittps : //dec . ligo.org/cgi-bin/DocDE/ShowDocument ?docid=2974| 
[37] S. Hild, S. Chelkowski, and A. Freise (2008), arXiv:0810.0604. 
[38] L. Blanchet, G. Faye, B. R. Iyer, and S. Sinha, Class. Quantum. 

Grav. 25, 165003 (2008), arXiv:0802.I249. 
[39] L. Blanchet, G. Faye, B. R. Iyer, and B. Joguet, Phys. Rev. D 

65, 06150I(R) (2002), Erratum-ibid 71, 129902(E) (2005), gr- 

qc/0105099. 
[40] L. Blanchet, T Damour, G. Esposito-Farese, and B. R. Iyer, 

Phys. Rev. Lett. 93, 091101 (2004), gr-qc/04060I2. 
[41] K. G. Arun, B. R. Iyer, B. S. Sathyaprakash, S. Sinha, 

and C. Van Den Broeck, Phys. Rev. D 76, 104016 (2007), 

arXiv:0707.3920. 
[42] C. Van Den Broeck and A. S. Sengupta, Class. Quantum Grav. 

24, 1089 (2007), gr-qc/0610126. 
[43] L. S. Finn, Phys. Rev. D 46, 5236 (1992). 
[44] L. S. Finn and D. F Chernoff, Phys. Rev. D 47, 2198 (1993). 
[45] C. Helstrom, Statistical Theory of Signal Detection, vol. 9 of 

International Series of Monographs in Electronics and Instru- 
mentation (Pergamon Press, Oxford, U.K., New York, U.S.A., 

1968), 2nd ed. 
[46] L. A. Wainstein and V. D. Zubakov, Extraction of Signals from 

Noise (Prentice-Hall, Englewood Cliffs, 1962). 
[47] C. Van Den Broeck and A. Sengupta, Class. Quantum Grav. 24, 

155 (2007), gr-qc/0607092. 
[48] L. P Grishchuk, V. M. Lipunov, K. A. Postnov, M. E. 

Prokhorov, and B. S. Sathyaprakash, Usp. Fiz. Nauk 171, 3 

(2001), astro-ph/0008481. 
[49] J. Abadie et al. (LIGO Scientific Collaboration), Classical 

Quantum Gravity 27, 173001(2010), arXiv: 1003.2480. 
[50] I. Mandel, D. A. Brown, J. R. Gair, and M. C. Miller, Astrophys. 

J. 681, 1431 (2008), arXiv:0705.0285. 
[51] A. Sesana, J. Gair, I. Mandel, and A. Vecchio, Astrophys. J. 

698, LI29 (2009), arXiv:0903.4I77. 
[52] C. Cutler and E. E. Flanagan, Phys. Rev. D 49, 2658 (1994). 
[53] N. Yunes, K. G. Arun, E. Berti, and C. M. Will, Phys. Rev. D 



18 

80, 084001 (2009), arXiv:0906.0313. [56] T. Damour and A. Nagar, Phys. Rev. D 81, 084016 (2010), 

[54] P. Ajith et al., Phys. Rev. D 77, 104017 (2008), arXiv:091 1.5041. 

arXiv:0710.2335. 
[55] Y. Pan et al. Phys. Rev. D 81, 084041 (2010), arXiv:09 12.3466. 



